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Abstract 

This work relates the theoretical study of the dynamic behavior of a dual free-piston Stirling engine (DFPSE) coupled with an asyn¬ 
chronous linear alternator. This machine integrates one piston and two displacers placed in a symmetrical way compared to the piston to 
improve the stability of the machine. The paper presents an analytical study of the dynamic balance equations of a DFPSE. This model 
takes into account the non-linear dissipative effects of the fluid and the electromagnetic forces. The dynamic balance equations of the 
machine are solved by means of linearized pressure in the time domain especially. The objective is to evaluate the thermo-mechanical 
conditions for stable operation of the engine. The developed model may be used to simulate the dynamic behaviour of a built engine. 
The DFPSE produces a mechanical power of 1 kW and it has a design operating point of 1.4 MPa corresponding to the frequency about 
22 Hz. Helium is the working fluid. This machine is designed to be used as a micro combined heat and power (gCHP) system for com¬ 
bined generation of electricity and heat. 

© 2006 Elsevier Ltd. All rights reserved. 
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1. Introduction 

Currently, there is a strong and media pressure for 
development systems of energy production and transfor¬ 
mation which are powerful with a durable policy of devel¬ 
opment [1]. The growing worldwide demand for less 
polluting form of energy has led to a renewed interest in 
the use of cogeneration technologies [2]. The concept of 
micro combined heat and power (pCHP) or micro cogene¬ 
ration has been known for long time. Cogeneration systems 
have the ability to produce both useful thermal energy and 
electricity from a single source of fuel such as oil or natural 
gas. 

There is an increasing interest in the use of Stirling engine 
based cogeneration systems for residential and commercial 
because of their prospect for high efficiency, good perfor¬ 
mance at partial load, fuel flexibility, low emission level, 
low vibration and noise level [3]. The Stirling engine is an 
external combustion reciprocating engine developed by 
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Robert Stirling in 1816. It uses an external source of energy 
to heat the gas located inside a cylinder (chamber). This gas, 
under pressure, expands when heated, driving a piston to 
perform work. The expanded gas volume, having released 
much of its energy, is cooled and compressed before the next 
heating cycle. Regenerator is the principal element of the 
Stirling engine. It is a porous environment which must 
absorb and reject heat during compression and expand gas 
volume. Stirling cycle is more near Carnot cycle when this 
element is effective [4,5]. Stirling engines are classified 
according to their arrangement [6]: the Alpha , Beta and 
Gamma. The Alpha configurations have two pistons in sep¬ 
arate cylinders connected in series by a regenerator, a heater 
and a cooler. Both the Beta and the Gamma configurations 
use the displacement pistons arrangement, but the Beta 
arrangement has the piston and the displacer in the same 
cylinder while the Gamma arrangement uses different cylin¬ 
ders. For this study, we considered a Beta configuration 
(Fig. 1). This paper presents in the first part the mechanical 
characteristics of a dual free piston Stirling engine (DFPSE). 
Then the mathematical model is considered stating assump¬ 
tions relative to the dynamical forces balance applied to the 
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Nomenclature 



A 

cross-sectional area for moving element, m" 

<5 

coefficient of real part of the polynomial equa- 

A' 

internal cross-sectional area, m 


tion 

Cud 

displacer gas spring damping per mass unit, 

1 

coefficient of imaginary part of the polynomial 


N s m 1 


equation 

Cup 

piston gas spring damping per mass unit, 

A P 

pressure drop, Pa 


N s m 1 

£ 

coefficient of imaginary part of the polynomial 

c d 

displacer load damping per mass unit, N s m -1 


equation 

c P 

piston load damping per mass unit, N s m -1 

y 

isentropic coefficient 

d 

diameter, mm 

7 

coefficient of real part of the polynomial equa- 

D 

damping coefficient per mass unit, s -1 


tion 

[D] 

matrix of damping coefficient 

•A 

calculation parameter 

E 

calculation parameter 


piston-displacer phase angle, ° 

f 

frequency, Hz 

CO 

angular velocity, rad s -1 

/cut 

cut-off frequency, Hz 

z 

coefficient of imaginary part of the polynomial 

F 

external force per mass unit, m s -2 


equation 

[F] 

matrix of external forces 



j 

imaginary number, V—l 

Subscripts 

K 

stiffness coefficient per mass unit, s -2 

alt 

alternator 

[K] 

matrix of stiffness coefficient 

b 

bounce space 

m 

mass, kg 

c 

compression space 

P 

pressure, bar 

ch 

charge 

Q 

heat transfer, W 

cool 

cooling 

Q(<o) 

polynomial equation 

d 

displacer 

Q 

power, W 

e 

expansion space 

r\, r 2 

displacer-piston stroke ratio, Yj Y p 

g 

gas 

r 

working gas (helium) constant, J kg -1 K -1 

h 

hot 

S 

calculation parameter 

H 

gas spring hysteresis losses 

t 

time, s 

Im 

imaginary part 

E 

temperature, K 

i 

initial condition 

V 

volume, m 

k 

cold 

y 

displacement, m 

mec 

mechanical 

kde 

expansion clearance space, mm 

P 

piston 

kpc 

compression clearance space, mm 

Re 

real part 



reg 

regenerator 

Greek symbols 

th 

total heat power 

a 

coefficient of real part of the polynomial equa- 

0 

initial condition 


tion 

kde 

expansion clearance space, m 

p 

coefficient of real part of the polynomial equa- 

Jpc 

compression clearance space, m 


tion 




mechanical elements and their thermofluidic equivalent 
elements. Optimized geometrical (volumes, masses) and 
thermophysical (pressures, temperatures) parameters are 
calculated to obtain a stable running of the machine 
designed to produce a mechanical power of 1 kW. 

2. Dual free piston Stirling engine (DFPSE) analysis 

2.1. Description of the DFPSE 

Stirling engines are beginning to stage a comeback to the 
market since the development of the modern free piston 
Stirling engine [7]. Free piston Stirling engine was invented 


by Beale in 1964 [8]. The free piston engine technology 
based on the Beta configuration was developed to alleviate 
the technical barrier posed by leakage problems. Free pis¬ 
ton Stirling engines are expected to eliminate mechanical 
contact, friction and wear, and provide tight sealing of 
the casing, thus requiring no mechanical maintenance dur¬ 
ing an operating lifetime of about 10 years [8,9]. Free pis¬ 
ton Stirling engine use working gas like gas spring in 
order to give an adequate movement to the various ele¬ 
ments of the engine (piston and displacer). 

A double free piston Stirling engine (DFPSE) arrange¬ 
ment is shown in Fig. 2. The engine acts as a vibrator gen¬ 
erator in which gas is compressed and expanded by 
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Fig. 1. Beta Stirling engine. 
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Fig. 2. DFPSE system. 


oscillating system of masses and gas springs. The engine 
appears as a close gas filled vessel with three moving parts: 
a piston and two displacers. The engine presents two heads 
of which each one is heated by an external heat sources. 
The two gas spaces enclosed by the displacers and the 
heads of the vessel are the expansion spaces which temper¬ 
ature can reach 400-700 °C. The working gas (helium, 
hydrogen, air) acts as the lubricant. The engine has a com¬ 
pression piston with double effect (the dual piston) and 
two displacers laid out in a symmetrical way compared 


to the dual one. The dual piston moving acts directly on 
each adjacent compression space. The two gas spaces 
enclosed by the displacers and the dual piston faces are 
water-cooled at temperature of 40-90 °C. The gradient 
pressures between hot and cold spaces strongly influence 
the dynamic behavior of the machine. The power strike 
is provided by the force enacted by the expanded gas in 
the two hot spaces on the displacers which are coupled 
to the dual power piston bys gas springs only. Electricity 
is generated by the interaction between piston and a linear 
alternator. This interaction causes an electromagnetic 
force with viscous effect which depends on the piston 
velocity. The engine works with a rating of natural oscilla¬ 
tions. In the DFPSE, there are no mechanical linkages 
coupling the pistons or displacers. The motions of those 
components follow working gas pressure variations. The 
free piston with the attached linear alternator can be 
tightly sealed to prevent the leakage of the working gas 
for a substantial period of time. The buffer volume (spring 
gas) present under the compression piston makes it possi¬ 
ble to place all the electric device of electricity generation. 
Hysteresis losses of cycle are reduced and we can consider 
a system which functions in trigeneration (electricity, heat, 
cold). 

The major advantages of free piston engines include 
input and output versatility, quiet operation, zero wear, 
zero maintenance, long life, ease of interfacing with the 
electric grid, continuous power operation and potential 
for high efficiency. The major disadvantage is the difficulty 
to stabilize the movements of the different elements like 
dual piston and displacers. The goal of this work is to find 
the mechanical and thermodynamical parameters which 
allow a correct operation of the DFPSE designed to pro¬ 
duce a mechanical power of 1 kW. 


2.2. Generalized analysis 

From a mechanical mass-spring analogy (Figs. 2 and 3), 
we performed a dynamic study of a free piston Stirling 
engine. Various gas volume of the engine behave like 
springs and shock absorbers. Compressive forces action 
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Fig. 3. Mechanical analogy. 


on piston faces and displacer faces are integrated into 
the forces balance. Fig. 3 shows the dynamic behaviour 
of a free piston Stirling engine with connection spring 
stiffness K and damping coefficient D between piston 
mass m p and displacer masses m d \ and m d2 . The external 
mass of the engine casing m c is higher than the other 
moving masses (piston and displacers) and the effects will 
be neglected in the mechanical forces balance. Then, the 
DFPSE may be described by a viscous damped three 
degree of freedom spring-mass system under external 
forces [10,4]. The dynamic equations may also be written 
as the sum of the spring factor matrix [K] and the damping 
factor matrix [D] [10]. The effects of the gravitational field 
have also been neglected because they are quite small com¬ 
pared to the effects of the pressures involved in the machine 

[H]- 

The system of coupled differential equations describes 
the dynamic behavior of the DFPSE (Fig. 3). Model out¬ 
puts available are evolution pressure, displacer and piston 
displacement, the electric charge, the thermal and electrical 
power. The motion of the system is completely described 
by the positions of the respective piston and displacers 
at any time t from the respective equilibrium positions. 
The dynamic behavior of a DFPSE is presented in the 
matrix form (Eq. (1)), a matrix of the stiffness coefficients 
(Eq. (2)) and a matrix of the damping coefficients (Eq. 
(3)). These matrices are written per unit of moving element 
mass of the machine. 


Newton’s second law yields: 


'h ' 

m 

= M 

£ 

+ [D] 

S’ ^• l 

+ [F«] 

'O 

_1 


_y&2 _ 


<N 





K pp 

^dlp 

^d2p 

T*pp 

Anp 

T>d2p 


^pdl 

^dldl 

^d2dl 

A>dl 

Alldl 

Ai2dl 


^pd2 

^dld2 

^d2d2 

Ppdl 

Alld2 

Al2d2 



F v (t) 0 0 

0 F dl (t) 0 
0 0 F d2 (t) 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


[F(£)] represents the matrix of possible external forces (i.e. 
piloting). In this study neither the piston, nor the displacers 
are controlled and F v (t) = 0, F d \(t) = 0 and F d2 (t) = 0. The 
piston is submitted to an electromagnetic force correspond¬ 
ing to the linear alternator function only. The pressure gra¬ 
dients existing between the different moving elements 
(piston and displacers) strongly influence the dynamic 
behavior of the machine. Each applied force has a compo¬ 
nent normal to the displacement surface pointing up. The 
variation of the volume located in the hot part (expansion 
space) is related to the displacement of the displacer while 
that located in the cold part (compression space) will be re¬ 
lated to the displacement of the piston. 

The differential model represents the damped dynamic 
system with three degrees of freedom (Figs. 2 and 3) [8,12]. 


l_ 

= H 

£ v 1 

+ [D] 

S' ^• l 

"a 

_i 

_j>d2_ 


_y di _ 


_jd2_ 


The equations of motion of the different elements are for 
the piston: 

OTpj)p = A p (P cX - P c2 ) - A’ v {P hi - P hl ) 

- (C palt + C H p)^p (5) 


the displacer 1: 

iB d \'y d \ — A d (P Q \ P d\) A r (P c \ P di) 

- (Qlalt + f^Hd 1 Ifdl (6) 


the displacer 2: 

iB d2 y d2 — A d (P c 2 P qi) T r (Pd 2 Pci) 

~ (Ol2alt + C H d2)jd2 (2) 


The dynamic equations take into account of the interac¬ 
tion between the alternator and the pistons (the electric 
charge causes a braking force). 

The DFPSE is divided into two parts (Fig. 3). The upper 
half-machine contains the expansion space ( K\,D \), the dis¬ 
placer {m d i) and the compression space ( K 2 ,D 2 ). The lower 




















































806 


J. Boucher et al. / Applied Thermal Engineering 27 (2007) 802-811 


half-machine integrates the compression space (K 3 ,D 3 ), the 
displacer (m d2 ) and the expansion space ( K 4 ,D 4 ). The dual 
piston m p separates these two parts. In order to approxi¬ 
mate fairly good the reality, the isothermal Schmidt analy¬ 
sis is applied in order to linearize the pressure in the 
different working volumes of the machine (Figs. 2 and 3). 
If the total gas mass contained in the machine is known, 
the principle of conservation makes it possible within the 
framework of the assumption of isotherm of gas volumes 
to determine the pressures P and Pp The pressure P corre¬ 
sponds to the mean pressure in the gas and P{ the initial 
pressure just before oscillation of the piston and displacers 
start (y d = Tp = 0)- These pressures appear in their general 
formulations: 




+ 




Fh A d J 7 de 

T k T k 


( 8 ) 

(9) 

( 10 ) 


The mass of gas m, the initial pressure Pi before oscillation 
start and the two parameters \jj and S are written for each 
part of the machine. P c h is the pressure of the gas charged 
in the engine at ambient temperature before oscillation 
start. This pressure integrates the S 0 parameter which cor¬ 
responds to the initial condition before start when the tem¬ 
perature T k is the same in the entire volume of the machine. 

The mean gas temperature T r in the two regenerator of 
the machine (upper and lower half-machine) will be evalu¬ 
ated by the average logarithm of the temperatures of their 
corresponding hot (r h ) and cold (T k ) ends: 


T h ~ T k 

ln(7yTk) 



2.3. Linearization of the pressure in the machine 

We consider the upper and lower half-machine volumes 
(Fig. 3). The expressions of the pressures in theses different 
volumes are linearized in order to be used in the dynamic 
analysis (Eq. (8)). Referring to Fig. 2, the different volumes 
variations are as follows: 

• For the volume of the upper half-machine: 


V bl = (A' pl - A' rl )y p 1 (12) 

Vci = A p i(y p i +J P ci) _ (An ~ A,-i )v d | (13) 

r el =A d l(y d i + J'dcl ) (14) 

• For the volume of the lower half-machine: 

V b2 = (A' p2 + A' r2 )y p2 (15) 

Vd = A v2 (y pc2 - y p2 ) + (A d2 - A [2 )y d2 (16) 

V e2 =A d2 (y de2 -y d2 ) (17) 


(y de i and y de2 ) and (y pc i and y pc2 ) are the heights of the 
expansion and compression volumes 1 and 2 at initial 
condition before oscillations, respectively. We can con¬ 
sider the following conditions of symmetry concerning 
the geometry of the DFPSE: A pi = A p2 = A p and the 
displacements of the pistons: y p i = y p2 = y p . 


For the volume of the upper half-machine, the mean 
pressure is 



(18) 

(19) 


with 


S^ = 


^pjpcl , Fkl , Erl ( Fhi ( ^dl^del 


= 


T ki 
ApT p 
T^Si 


+ 


+ — + 


+ 


Tki T r i 
A r l ~ A d l 
T k \S\ 


+ 


+ 


T hi 

Adi 

TmSi 


lTdi 


( 20 ) 

( 21 ) 


The mean temperature of the regenerator of the upper half- 
machine is 


T _ F hl - T kl 
rl ~ln(T hl /T k l) 

Then, after a binomial expansion of [1 + i/q] _1 , neglecting 
the second-order terms, the linearized pressure becomes: 



P i = P a 


Apy v 

T kl S 


Ad i — A 


rl 


T k i^i 


A d i \ 
TmSi) 


Tdi 



Pci ~ P\ 

- A 

P eX « Pi + A P x = P\ + ^fy p + 


(24) 

(25) 

(26) 



Pdl 

Adi 


Tdi 



For the lower half-machine, the mean pressure is 



with 




Apy p c2 

T kl 

Apy v 

T k2 S 2 


+ 

+ 


F k2 Fr 2 Fh 2 A d2 y de2 

rji 1 rri 1 rji 1 rri 

I k2 ‘ r2 ‘ h2 ‘ h2 


A r 2 — A d2 A d 2 

+ \y d2 


T k2 S 2 


T h2 S 2 


(28) 

(29) 


(30) 

(31) 


The mean temperature of the regenerator of the lower half- 
machine is 


T 1,2 — Tki 

In (T h2 /T k2 ) 
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Then, after a binomial expansion of [1 + \j/ 2 ] \ neglecting 
the second-order terms, the linearized pressure becomes: 



Ad2 ~ A 


r2 


T k2 S 2 


^d2 \ 

T\7s~ 2 J yd2 


Pb2 = Pli f 1 - 
P C 2 ~ P2 
Pdi = Pn U - 


K 2 

V b2l 



yArl 

V d 2i 


}’d2 


(33) 

(34) 

(35) 

(36) 


P 


e2 ~ 


P 2 + AP 2 =P 2 - 


A 


d2 


c d2 . 

T-Td2 

^d2 


E 


p 2 


A 


y 


d2 


P 


Ed2 
Ad2 


y d 2 



matrix [D] is composed by elements which generate equiv¬ 
alent viscous losses. C pa i t corresponds to a loss of energy 
caused by a viscous dissipation of electromagnetic origin 
in the electromagnetic coil. C Hp , C H di and C Hd 2 are hyster- 
etic loss due to cyclic flows between the different elements 
(piston and displacers). These losses depend on the mass 
flow generated by the piston movements, the frequency 
and the magnitude of the pressure [13-17]. In this work 
the two displacers 1 and 2 are not controlled. Then, the 
corresponding damping coefficients do not exist and 

f^dlalt 0 and C<d2alt 0. 

2.5. Analysis of the harmonic displacements of the pistons 
and displacers 


2.4. Stiffness [K] matrix and damping [D] matrix 

The dynamic behavior of a FPSE is presented in the 
matrix form with a matrix of the stiffness coefficients as 
well as a matrix of the coefficient of dissipations. We con¬ 
sider the dynamic equations for a thermomechanical sys¬ 
tem with three degrees of freedom. Eqs. (5)-(7) may be 


presented in the matrix form: 

Tp = ^pp3p "E ^pdiJdi "E K-pd2ydi ~ ^ppTp (38) 

Jdl —^dlpTp + ^dldlJdl "E T^dipjfp T^dldljfdl (39) 

y &2 = Kd 2 V y v + Kd 2 d 2 yd 2 + ^d 2 p jp — T*d 2 d 2 jd 2 (40) 

with 

K vv K v d\ K V d2 

[K] = ^di P *didi 0 (41) 

_^d2p 0 Kd2d2 _ 

~ D vv 0 0 

D] = An p Audi 0 (42) 

_Al2p 0 Dd2d2 _ 


Different assumptions relate to the pressure within volumes 
and the initial pressure (Fig. 2). 

The initial mean pressures within the upper and lower 
volumes are equal: 

Pn = P 2 i = Pi (43) 

The surface areas of the piston and the two displacers are 
equal: 

Ad 1 = Ad2 = Ad (44) 

The internal surface areas of the two displacers are equal: 
A r 1 = A r 2 = A r (45) 

The displacements of the piston and displacers are equal: 

Jpcl = y p cl = Jpc (46) 

3de I 3\!e2 3de (47) 

The linearized coefficients of the stiffness and damping 
matrices are given by the Table 1. The damping load 


We shall consider the piston and displacers 1 and 2 to be 
harmonic. Complex notation is used to describe the 
response of the components: 

y p =y p ^\ y dl =y di e^ l \ 

y d 2 = y d 2 e i( “ + " 2) ( 48 a) 

The time derivatives are: 

>p = j e °> ydi = j 0 )y d , e> 1 P ,+<p1 >, 

y<s 2 = j® 3 ; a 2 e i( “' +<?2) ( 48 b) 

j) p = -co 2 y p e (j m t \ j> d i = —a> 2 y dl 

y d 2 = -(o 2 y d 2 e (t0 ' l+(p2) ( 48 c) 

Substituting Eqs. ( 48 a)-( 48 c) into Eqs. ( 38 )-( 40 ) yields: 
K pp + nKpii e ^ 1 + r 2 K pd2 e ^ 2 + jcoZ) pp + co 2 = 0 ( 49 ) 

T^dip + r l^diai + j^Aiip + j ojr l-Ddidi e ^ 1 + <x> 2 r\ = 0 

( 50 ) 

Al2p + r 2 K d 2 d 2 d ,p2 + )<x>D dlv + ]wr 2 D d 2 d 2 P p2 + orr 2 d ,pl = 0 

( 51 ) 


We define: 


• the ratio of the maximum displacement r\ between the 
piston and the displacer 1: 


r\ = 


Tdi 


max 


y 


(52a) 


pmax 


r 1 


{CO (X(jip “I - ^dlp^dldl ) “^-^dldl^dlp) CD~ {CD"DE^\^\D X^dlp-^dldl) 


(co 2 +Xdidi) + (co 2 D d idi) 


(52b) 

• the ratio of the maximum displacements r 2 between the 
piston and the displacer 2: 


r2 = 


r 2 


Td2 


max 


y 


(53a) 


pmax 


(c0 2 (K d2p + Ai2pAi2d2) +^d2d2^d2p) + CO 2 ( 03 2 £>d 2 p +^d2d2^d2p ~ K& 2 pD fo&lY 

(CO 2 + K& 2 & 2 .Y + (m 2 Z)^ 2 d 2 ) 


(53b) 
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Table 1 

Elements of the stiffness and damping matrices 

Elements of the stiffness matrix Elements of the damping matrix 


K 

K 

K 

K 

K 

K 

K 


1 


1 \A: 


pp 


+ — 

8 2 T k S\T k J m v 


2 y A l 




Pi 


p V b' 


_ A p f A d \ A r \ A d i 

pdl i 


m 


p 


SiT, 


S\T 


hi 


_ Ap f A d2 — A r2 A d2 

pd2 i 


m. 


S 2 T k 


S 2 T 


2 1 h2. 


dip 


A p A r \ D 1 | E p i 

■ O — — r 


m d i 


T k S { ntdi 


d2p 


A p A r2 


Pi 


1 


+ 


'p2 


A T 1 

dim =—P, 

m d i 


A r 2 

_ rz D 

d2d2 — -“i 

m d 2 


fU d2 m d2 

Aril — A rl T 


T k S\ 

Th\Si 

Vm) 

i + 

m d i 

d2 — A r 2 

A d2 

yA r2 \ 

1 + 

E d2 


T A 2 T k 2 S 2 V d02 


m d2 


Dpp — -(Cpait + C Hp ) 

m n 


T> p di — 0 


^pd2 — 0 


f^dlp 


c 


pi 


m d i 


D d2 p 


C 


p2 


W d 2 


Alldl — -(Oil _ E*Hdl — Ed l a it) 

m d 1 


Al2d2 — -(C d2 — E H d2 — Ed2alt) 

m d2 


the phase cp\ between the piston and the displacer 1: 


r i e 


(j<pi) _ _ 


(j) l = tan 


-l 


Ailp A~ jfuAllp 
CO 2 + Alldl + j&Alldl 

Co(n) Ailp T~ AlldlAllp -Allp-Alldl) 
E^CAllp "h AilpAildl) A~ -Andl-Allp 


(54a) 


(54b) 


the phase cp 2 between the piston and the displacer 2: 


r 2 e 


(j>2) _ _ 


02 = tan 


-l 


^d2p T jmAl2p 
E) 2 + Al2d2 + j<T>Al2d2 

Co(cO F) d2 p + -Ai2d2'Al2p -Al2pAi2d2) 

E^ 2 (^d2p H” ^d2p-^d2d2 ) H" ^d2d2^d2p 


(55a) 


(55b) 


The solution of the above system of differential Eqs. (49)- 
(51) is based on the roots of the polynomial equation 
Q(co): 

Q(co) = — (co 2 A- K pp A~ jcoD pp )(co 2 + Audi + j<r>Aiidi) 

X (co 2 + Al 2 d 2 + j&Al 2 d 2 ) + K pd \ (Allp A~ }(oD d \ p ) 

X (co 2 T d 2 d 2 A- )(OD& 2 & 2 ) + ^pd 2 (Al 2 p + j&Al 2 p) 

X (CO 2 + Alldl + j E0Z>d 1 d 1) = 0 (56) 


QH = 0 Re (co) + j Q lm (co) (57) 

g Re (co) = (co 6 — /ico 4 + <5co 2 — a) (58) 

0imM = (x™ 5 - ^ - &) (59) 


Sixth power of pulsation shows many oscillations modes 
are possible [10]. We study the auto excited periodic modes. 
For the free oscillations characterized by no external forces 
applied on the piston and on the two displacers (F p (t) = 0, 
F d \(t) = 0 and F d2 (t) = 0), the real part of the characteristic 
polynomial g Re (co) is always of an order of power higher 
(6th degree) than the imaginary part of the characteristic 
polynomial Qi m (co) (5th degree). Steady oscillations appear 
when the characteristic polynomial has two imaginary 
roots and four roots with negative real parts [10,12]. Then, 
for a symmetrical machine in steady oscillations Qi^co) = 
0 and the characteristic polynomial Q(co) satisfy the follow¬ 
ing constraints [12-14]: 

a Audi T ^pp ^ 0 

P = A^P Alldl T AlldlA>P ^ ^ 
y Aiidi-App F) pp K d \ d i T 2An P Adi ^ ^ 

3 — ^pp^Alldl ^^dlp^pdl ^ 0 

= o 

a y 


(56a) 

(56b) 

(56c) 

(56d) 

(56e) 


This characteristics equation determines the 6 independent 
roots for the damped vibration problem and it can be put 
in the general complex form Q(co) defined by the real part 
2 Re (co) and the imaginary part Qi^co): 


3. Results 

A theoretical parametric study of a Stirling engine with 
a dual piston and two displacers is conducted to obtain the 
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frequencies and the operating capacities. The mechanical 
power can reach 1 kW for a frequency of 20-35 Hz. During 
dimensioning phase, internal dynamics of engine is actually 
unknown since no system mechanical imposes neither dis¬ 
placement of pistons, no phase lag between piston and dis¬ 
placer (on the contrary traditional Stirling engines whose 
rotational movement is recovered on the linkage) [18]. 

The matrix coefficients depend on the geometry of the 
engine. A linearized dynamical model, with stiffnesses 
matrix [K] are fixed by the engine geometry and the average 
pressure gas. The three parameters of the damping load 
matrix Z> pp , Z> dldl and D d2 d 2 depend on the constant elec¬ 
tric charges and thermal dissipation coefficients in gas vol¬ 
ume. E d i, E d2 and E p are corrective terms which multiply 
the pressure losses depending on the velocity and the dis¬ 
placement. The DFPSE characteristics are given in the 
Table 2. 

Figs. 4-10 show the behavior of the FDSPE at stable 
operation with the different characteristics defined in Table 
2 (Figs. 4 and 5). The charge pressure and the hot source 
temperate determine the mechanical power produced, the 
frequency of operation as well as the stroke of the dis¬ 
placer. The amplitude of the stroke and the phase between 
displacers and piston vary with the pressure. The strokes of 
the piston and displacers are 2 cm and 1 cm, respectively. 
The mean pressure in the machine is 1.4 MPa for a 
dynamic amplitude of the pressure of 0.1 MPa. 

The curves of the Figs. 6-10 present the evolution of the 
mechanical power produced by the piston, the thermal 
power, the stroke of the piston, the amplitude ratio of 
behavior displacers and piston and finally the phase angle 
between the pistons and the displacers. The stable opera¬ 
tions appear between 20 Hz and 35 Hz. For frequencies 
lower than 20 Hz, the operation of the engine is unstable 
and the amplitudes of piston and displacers grow strongly 
with unrealistic values. For frequencies higher than 35 Hz, 


Table 2 

DFPSE characteristics for the simulation 


Dual free Stirling piston engine data 


General 

Working fluid 

Helium 

Gas constant 

r = 2080 J kg -1 K -1 

Specific heat 

c p = 5183 J kg -1 KT 1 

Isentropic coefficient 

y — 1.667 

Dynamic viscosity 

/i = 1.8 x 10~ 5 PI 

Charged pressure 

P ch =10 bar 

Mean cold space temperature 

F k = 320 K 

Mean hot space temperature 

71= 815 K 

Geometric 

Piston diameter 

<7 P = 120 mm 

Internal piston diameter 

d p i =116 mm 

Diameter of the displacer rod 

d r — 35 mm 

Displacer diameter 

d d dp 

Piston mass 

m p — 6.2 kg 

Displacer mass 

m d = 0.426 kg 

Height of the compression space at rest 

y pc =150 mm 

Height of the expansion space at rest 

Vde =100 mm 


Stroke (m) /Piston 



Fig. 4. Piston stroke for a frequency of 33 Hz. 

it may be seen that it is impossible to extract power from 
the engine. For these frequencies the power losses by fric¬ 
tion become very important and the different amplitudes 
of piston and displacers are damped until the mechanical 
power tends to zero. Fig. 6 shows the mechanical power 
of the engine versus frequency. With the increase of fre¬ 
quency, the mechanical power firstly increases to the max¬ 
imum and then decreases gradually until zero power for the 
frequency of / cut = 37 Hz which correspond to a cut-off 
mechanical power frequency. The maximum mechanical 
power is Q mec = 1050 W and the corresponding frequency 
is /= 22.5 Hz. Fig. 7 presents the thermal power of the 
machine versus frequency. This power corresponds to the 
input thermal power by an external source at one side of 
the machine only. The total input heat power of the 
FDSPE is the double. The maximum power decreases with 
the frequency and presents a maximum value Q h = 1750 W 
for the frequency / =20 Hz. The total heat power is 
2th = 2gh = 3500 W. The cut-off frequency occurs for 
/cut = 37 Hz. Fig. 8 shows the behaviour of the cooling 
power versus frequency. This power decreases with the fre¬ 
quency in the same order than the mechanical and thermal 
powers. The maximal cooling power is Q coo \ =2500 W for 
the frequency / =20 Hz. The cut-off frequency occurs for 
/cut = 37 Hz. Figs. 6-8 show the conservation of the power 
in the system and: 

2th = 22h = Scool + 2mec (37) 

Fig. 9 presents the displacement ratio versus frequency 
which is bell-shaped in the range of stable operation 20- 
37 Hz. The ratio of maximum displacement is r\ = y d 2 max/ 
kpmax = 0.61 for the frequency / = 24.5 Hz. This ratio of 
displacement does not appear for the same frequency than 
for the mechanical and thermal powers. This ratio of dis¬ 
placement varies from 0.52 to 0.61. The amplitude of dis¬ 
placement of the displacer is increasingly smaller than 
that of the piston. Fig. 10 indicates the relation between 
the phase angle between piston and displacers and the fre¬ 
quency. The movements of the displacers are supposed to 
be identical. The machine is built in a symmetrical way 
and we suppose cp\ = (p2 = <P- The phase angle decreases 
with the frequency and the optimum value (p = —60° (or 
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Fig. 5. Pressure amplitude for a frequency of 33 Hz. 


Q mec (w) 



Fig. 6. Mechanical power versus frequency. 


Qcool (W) 



Q,h W 



cp = 300°) is obtained for 21 Hz. The negative value of the 
phase angle indicates piston motion lagging behind those 
of the displacers. For the stable operation of the machine, 
the piston-displacers phase angles lie in the fourth quad¬ 



rant (180° < cp\^ 2 < 360°) depending of the size of the 
denominators in Eqs. (55a) and (55b). For a phase angle 
of 0°, it corresponds to a cut-off frequency of 35 Hz for 
which the DFPSE do not produce work. 
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<p(°) 



/(Hz) 


Fig. 10. Phase cp between piston and displacers versus frequency 
Ol = (P 2 = (P)- 


4. Conclusions 

The dynamic behaviour of a dual free-piston Stirling 
engine (DFPSE) coupled with an asynchronous linear alter¬ 
nator has been performed with a mathematical analysis 
including complex damping and stiffness coefficients. The 
dynamic balance equations of the machine are solved by 
means of linearized pressure in the time domain especially. 
The objective is to obtain a stable operation. This machine 
integrates one piston and two displacers placed in a sym¬ 
metrical way compared to the piston to improve the 
mechanical stability of the machine. We studied the design 
of a DFPSE to produce a mechanical power of 1 kW. The 
stability of a DFPSE is obtained taking into account vari¬ 
ous efforts like the geometric, dynamic and thermodynamic 
variables. Thermo-mechanical simulation makes a first 
study useful for the predimensioning of the free piston Stir¬ 
ling engine. We thus obtain from values geometrical an 
operating range of frequency, recovered electric power, 
the amplitude ratio of displacement of the pistons, the vari¬ 
ations of the pressure in the machine. The stable operation 
was performed for a frequency in the range 20-35 Hz. The 
strokes of the piston and displacers are 2 cm and 1 cm, 
respectively. The mean pressure in the machine is P mean = 
1.4 MPa for a dynamic amplitude of the pressure of 
AP = 0.1 MPa. 
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